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I. Introduction 


In this project, research has focused upon studies for the heavy-ion-dominated plasma torus 
(0 + , 0 ++ , S+, S ++ , S +++ , S ++++ ) and its near plasma sheet extension produced by the ionization 
of neutral gases from Io. Specifically, we have undertaken explicit studies and calculations to 
describe the spatial nature of the Iogenic plasma source and have begun to study its outward 
rotation-driven transport pattern in Jupiter's magnetosphere in order to identify and understand 
the space plasma phenomena that shape it into the Io plasma torus. The approach has been to 
combine state-of-the-art calculations at AER for the Iogenic plasma source and numerical 
solutions of the nonlinear transport equations in a collaborative effort with Drs. T. W. Hill, R. A. 
Wolf and colleagues at Rice University using the Rice Convection Model for Jupiter (RCM-J). 
In this way, we have begun to address the complex plasma torus structures that emerge from the 
source and transport processes. 

The RCM-J's use of a fine two-dimensional grid and its many-fluid representations of the 
plasma make it capable of representing this complex region. The calculation of a 
continuous/spacetime-dependent Iogenic plasma source in three-dimensions and its proper 
interface with the RCM-J provide a critical and, up to now, missing element for this convection 
model. The Iogenic plasma source model is quite sophisticated and contains, for example, east- 
west and System III longitudinal asymmetries in the spacetime description of the plasma torus. 
The many different observational data that are carefully crafted together in determining the 
nature of the Iogenic plasma source have brought a required reality to this convection model that 
is necessary to unravel the structure and transport physics of the plasma torus. These numerical 
simulations are ultimately intended both to clarify the conditions for validity of various 
theoretical scenarios and to provide a useful tool for interpretation of plasma torus structures and 
in situ and remote data acquired by Galileo. 

II. Overview of Work Performed in the Project 

2.1 Iogenic Plasma Source 

The Iogenic plasma source can be divided into two regions. First, there is a very spatially 
confined and denser "Inner Region" below Io's exobase where pickup (i.e., both electron impact 
ionization and charge exchange) ions are created in the local atmosphere and also where 
Pedersen and Hall conductivities become important in the collisionally bound ionosphere. 
Second there is a spatially extended but less dense "Outer Region" above the exobase created by 
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neutral gases from Io in the bound corona and in the more distant escaping neutral clouds where 
pickup ions are also created. The solution of the Iogenic plasma source for the small Inner 
Region" is very complex and involves the coupling of the atmospheric structure to the 
electrodynamic interaction problem. The recent modeling of Saur et al. (1999) explicitly 
concentrated on this Inner Region with the primary purpose of solving the electrodynamic 
problem but with a fixed spherically symmetric atmosphere. The strength of the Iogenic plasma 
source for the Inner Region (effectively a point source on the broader circumplanetary scale) is 
thought to be comparable to and even likely larger than the Iogenic plasma source for the Outer 
Region. It is the volumetrically-concentrated Iogenic plasma source for the Inner Region that is 
now recognized (Smyth and Marconi 1998; Hill and Pontius 1998; Delamere et al. 2001) to be 
the dominant plasma source for creating the main density peak in the plasma torus, the so-called 
plasma ribbon, that is located at western elongation at -5.7 Jupiter radius. 

The primary research undertaken in this project has been to investigate the 
circumplanetary nature and the near-Io nature of the Iogenic plasma source for the Outer Region. 
The general nature of the Iogenic plasma source for the Outer Region was deduced earlier by 
Marconi and Smyth (1996). They noted that the nominal lifetimes (-1 to 20 hrs) for neutrals (H, 
O, S, SO, S0 2 ) in the plasma torus [due to electron impact ionization (or dissociation) and/or 
charge exchange reactions] are short compared to the characteristic time of -200 hrs for neutrals 
to form a complete 360° toroidal cloud about Jupiter at Io's orbit. This rapid destruction of the 
neutral clouds in the torus hence produces a spatially extended but rather circumplanetary 
asymmetrically Iogenic plasma source both upstream and downstream of Io that has a large peak 
at the satellite’s location. In this project, the Iogenic plasma source for the Outer Region has 
been calculated using the AER Neutral Cloud Model. The spatial nature of the Iogenic plasma 
source created by atomic oxygen and atomic sulfur species in the Outer Region has been 
calculated in three dimensions and investigated in detail both on a circumplanetary scale and on 
a near-Io scale for one location of Io in local time. These calculations are discussed in Section 
III and are presented in two companion papers: “Nature of the Iogenic Plasma Source in Jupiter s 
Magnetosphere I. Circumplanetary Distribution” (in appendix; Smyth and Marconi 2001a) and 
“Nature of the Iogenic Plasma Source in Jupiter’s Magnetosphere II. Near Io Distribution” (in 
preparation, Smyth and Marconi 2001b). 
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2.2 RCM-J Transport Calculations 

The original RCM-J model (Yang 1992; Yang et al. 1992, 1994; Pontius et al. 1998) has 
been improved during this project with the collaborative project at Rice University so as to 
include a space-time dependent Iogenic plasma source. This inclusion of the Iogenic plasma 
source has required a major revision and testing at Rice University of the earlier model code. In 
the improvement of this model, many unforeseen difficulties have been encountered. Although 
significant progress has been made, the advancement has been much more difficult and slower 
than initially expected. 

The spatial nature of the Iogenic plasma source developed in this project has been made 
available to Rice University for the purpose of being incorporated into the RCM-J transport 
calculations. For the stage of development that the RCM-J has achieved during this project, 
however, the three-dimensional circumplanetary and near-Io distributions developed at AER 
contain too much spatial realism for the present state of the model to incorporate, although 
progress in the past year has been moving forward at a faster pace. A summary of the new 
calculations incorporating various physical elements of the Iogenic plasma source is presented in 
Section IV. 

III. Model Calculations for the Iogenic Plasma Source in the Outer Region 
3.1 Introduction 

Efforts have been focused upon the investigation of the nature of the Iogenic plasma 
source for the Outer Source on two different spatial scales: (1) on a circumplanetary Jovicentric 
scale including all of Io's circular orbit (radius of ~5.9 Rj) around Jupiter and covering the radial 
interval from 4 to 10 Rj, and (2) on a near-Io scale for distances from the nominal exobase of 1 .4 
Io radii (Ri 0 ) radially outward including the Lagrange sphere (radius -5.85 R l0 ) to a distance of 
~±20 R Io [i.e., ~±0.5 Jupiter radii (R,)] from Io. These calculations have been undertaken using 
the AER Neutral Cloud Model. Efforts have been focused upon the Iogenic plasma source that 
are created only in the Outer Region by O and S atoms ejected from Io’s exobase. The 
calculations are restricted to O and S since they are known to exist and have exobase source 
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estimates that are more certain. In these calculations, the incomplete collisional cascade source 
velocity distribution that was deduced for Na at Io's exobase by Smyth and Combi (1997) has 
been adopted for the shape of the O and S source velocity distributions at Io's exobase. An 
isotropic source rate at Io’s exobase of 3.2 x 10" 7 atoms s 1 was initially adopted for O at an 
assumed exobase radius of 2600 km and was determined from energetic monoenergetic 
modeling efforts of Smyth (1992). The S source rate was initially assumed to be 1.6 x 10 27 
atoms s' 1 , half of the O source rate, as would be consistent for a parent molecule S0 2 . Later in 
the project, however, more realistic O and S source rates of 1.38 x 10" 8 s 1 and 0.69 x 1 0 28 s 1 , 
about ~4 times larger than the earlier rates, were adopted. These new rates were based upon the 
more recent modeling of O neutral clouds by Smyth and Marconi (2000) that adopted for O the 
incomplete collisional cascade source velocity distributions deduced for Na at Io's exobase by 
Smyth and Combi (1997). All calculations discussed below were undertaken for Io located at 
western elongation (i.e., a geocentric phase angle of 270°, which is assumed to be at the local 
time of the dusk ansa) and at an Io System III longitude of 230°, so that it is well north of the 
centrifugal equator symmetry plane of the plasma torus and has a magnetic off-set tilted dipole L 
value of 5.88. 

3.2 Summary of the Spatial Nature of the Iogenic Plasma Source: Circumplanetary Scale 

Calculations have been performed on a circumplanetary scale about Jupiter to study the 
asymmetric spatial longitudinal and radial asymmetric nature of the O and S neutral densities and 
the ion pickup rates that they create in the proximity of Io's orbit. For this large spatial scale, the 
O and S neutral densities have been integrated perpendicular to Io’s orbit plane to produce a two- 
dimensional distribution with a spatial resolution box size of (2 x 10 4 ) x (2 x 10 4 ) km. The 
instantaneous ion pickup rates have been integrated along the magnetic field lines to produce a 
two-dimensional distribution in the plasma torus centrifugal plane with an angular longitudinal 
resolution of 5° and a radial resolution of 0. 1 Rj in the radial interval from 4 to 10 Rj. 

Only a brief summary of the paper “Nature of the Iogenic Plasma Source in Jupiter’s 
Magnetosphere I. Circumplanetary Distribution” is given here, since the paper is included in the 
appendix. Two-dimensional distributions produced by integrating the three-dimensional 
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information along the magnetic field lines are presented for the pickup ion rates, the net-mass 
and total-mass loading rates, the mass per unit magnetic flux rate, the pickup conductivity, the 
radial pickup current, and the net-energy loading rate for the plasma torus. All of the two- 
distributions are highly peaked at Io’s location and hence highly asymmetric about Jupiter. The 
physical origin and the amount of asymmetry for these volumetric quantities are discussed. The 
spatially integrated rates for the net-mass and total-mass loading rates are 154 kg s 1 and 277 kg 
s 1 , respectively, with 71% of the rate confined to a very large peak located within ±20° longitude 
of Io’s orbital location. The spatially integrated rate for the mass per unit magnetic flux is 1 . 1 x 
10 8 kg s’ 1 T 1 . Maximum values for the pickup conductivity and the radial pickup current at Io s 
location are 2.5 x 10' 3 mho and 4 x 10 4 amps. The spatially integrated net-energy loading rate is 

1.03 x 10 12 watts and is dominated by the contributions from the peak centered on Io. 

The space-time variability of the Iogenic plasma source is discussed and is shown to 
depend upon the complete time history of the neutral clouds from which they are created. The 
Iogenic plasma source calculated here for the Outer Region is compared to various estimates for 
the Iogenic plasma source created in the combined Inner and Outer Regions. The mass source 
rate for the Inner Region is estimated to be -300-700 kg s' 1 , about 1 to 2.5 times larger than that 
calculated for the Outer Region (277 kg s '). Other possible sources of pickup ions for the Outer 
Region are also discussed. The Inner Region is therefore expected to be the dominant 
contributor of ions to the formation of the plasma torus ribbon, the maximum density peak in the 
radial structure of the plasma torus, and to be a comparable if not larger contributor to the energy 
budget of the plasma torus. Future improvements in calculating the Iogenic plasma source 
include the addition of near-Io plasma effects in the satellite flank and plasma wake and the more 
complete specification of source rates for neutrals at Io’s exobase. 

3.3 Spatial Nature of the Iogenic Plasma Source: Near-Io Scale 

Since the Iogenic plasma source is shown in Section 3.2 to be very highly peaked at Io's 
instantaneous position, we have undertaken near-Io scale calculations of much higher spatial 
resolution to investigate the behavior of the large spatial gradients near the satellite. Calculations 
presented below are for the initial neutral source rates noted in Section 3.1. Calculations in the 
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paper “Nature of the Iogenic Plasma Source in Jupiter’s Magnetosphere II. Near-Io Distribution’ 
will be similar to results presented here but will adopt the new and more accurate neutral source 
rates noted in Section 3. 1 . For the O and S isotropic source velocity distributions adopted above, 
the spatial nature of the Iogenic plasma source has been calculated in three dimensions and 
collected within spatial grid boxes contained within two different volumes centered around Io's 
instantaneous location. The linear dimension of the smaller volume about Io is ±1.05 x 10 4 km 
(or ±5.785 R Io ) and is composed of grid boxes with a finer spatial resolution of 600 km (or 0.331 
R l0 ). This volume is comparable to the volume of Io's Lagrange sphere, which has a radius of 
5.85 R Io . Along Io's orbit, this smaller volume is equivalent to a longitude angular interval of 
only ±1.425° about the satellite. The linear dimension of the larger volume about Io is ±4.375 x 
10 4 km (or ±24. 1 R l0 ) and is composed of grid boxes with a coarser spatial resolution of 2500 km 
(or 1.377 R Io ). This larger scale is important to investigate asymmetries that occur in the Iogenic 
plasma source because of spatial gradients in the plasma torus that exist near Io's orbit. Along 
Io's orbit, this linear dimension of this larger volume is equivalent to a longitude angular interval 
of ±1 1.889° about the satellite. 

3.3.1 Instantaneous Source Rates: Near-Io Scale 

For the atomic oxygen source, the model calculation of the instantaneous ion pickup rates 
for electron impact and charge exchange is shown, respectively, in Figures la and lb for a two- 
dimensional slice in a plane that contains Io and that is perpendicular to the Jupiter-Io outward 
radial coordinate. In Figure 1 at distances from Io larger than the Lagrange sphere, both the 
electron impact and charge exchange rates exhibit a north-south asymmetry, as expected. The 
larger rates below Io are created by the larger electron density in the plasma torus experienced by 
neutrals that are south of Io and hence closer to the plasma torus centrifugal equator. A 
longitudinal asymmetry in the downstream direction would also be expected if a description for 
the Io plasma wake had been included in the model calculation. 

Because of the steep radial gradient in the plasma torus near Io's orbit, the Iogenic plasma 
source will also exhibit a radial asymmetry. This radial asymmetry is exhibited in the profile 
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LONGITUDINAL DISTANCE FROM IO (R| 0 ) LONGITUDINAL DISTANCE FROM IO (R|o) 

Figure 1. Iogenic Plasma Source Rates for Atomic Oxygen. In a plane through Io and 
perpendicular to the Jupiter-Io outward radial coordinate, the two-dimensional spatial structure of 
the ion pickup rate (ions cm' 3 s' 1 ) that is produced by atomic oxygen is shown in (a) for electron 
impact ionization and in (b) for charge exchange reactions. Io's location and size are shown by 
the shaded disk. In the model calculation, an incomplete collisional cascade source of 3.2 x 10“ 
oxygen atoms s 1 was assumed at Io's exobase. The ion pickup rates were collected in coarser 
grid boxes with a spatial resolution of 2500 km. 

calculated along the Jovicentric radial coordinate through Io (i.e. perpendicular to the plane in 
Figure 1) and is shown, respectively, for both atomic oxygen and sulfur in Figures 2a and 2b. 
Beyond the Lagrange sphere distance, all the profiles inside of Io’s orbit can be seen to decrease 
more rapidly than the profiles outside of Io's orbit because of the steep gradient in the electron 
density and temperature in the plasma torus near Io's orbit. For atomic oxygen in Figures la, lb, 
and 2a, the charge exchange rate can be seen to be consistently larger (~3 or greater) than the 
electron impact ionization rate (El). In contrast, for atomic sulfur in Figure 2b the electron 
impact ionization rate (El) can be seen to be larger (—1.6) than the charge exchange rate outside 
of Io's orbit and comparable to the charge exchange rate inside of Io’s orbit where the electron 
temperature decreases. This difference in behavior occurs primarily because of the larger 
ionization potential for atomic oxygen (13.62 eV) as compared to atomic sulfur (10.36 eV). 
Within Io’s Lagrange sphere in Figures 1 and 2, the Iogenic plasma source rates are essentially 
symmetric about Io. For the finer grid boxes with a resolution of 600 km, an improved model 
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Figure 2. Iogenic Plasma Source Rates for Atomic Oxygen and Sulfur. Calculated radial 
profiles for the ion pickup rates (ions cm 3 s' 1 ) near Io for electron impact ionization (El) and for 
charge exchange (CE) are shown in (a) for the atomic oxygen source and in (b) for the atomic 
sulfur source. The location of Io is shown by the shaded area and that of the Lagrange sphere by 
the vertical dashed lines. An incomplete collisional cascade source at Io's exobase of 3.2 x IO 27 
atoms s ' and 1.6 x IO 27 atoms s' 1 was assumed, respectively, for oxygen and sulfur. The pickup 
rates were collected in the coarser grid boxes with a spatial resolution of 2500 km. 


calculation of this spatial gradient within Io's Lagrange sphere is presented in Figures 3a and 3b 
for the same calculation shown in Figures 2a and 2b. From Figure 3, the maximum 
instantaneous CE rate near Io's nominal exobase radius of 2600 km is for atomic oxygen about 
0.7 cm' 3 s' 1 and 2.1 cm' 3 s' 1 (for a ratio of 1 to 3) and for atomic sulfur about 1.8 cm' 3 s' 1 and 1.1 
cm' 3 s ' (for a ratio of 1.6 to 1). The combined oxygen and sulfur maximum pickup El and CE 
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Figure 3. Iogenic Plasma Source Rates for Atomic Oxygen and Sulfur. Same figure caption 
as for Fig. 2 except that there are no vertical lines for the Lagrange sphere and the pickup rates 
were collected in the finer grid boxes with a spatial resolution of 600 km. 

rates near the exobase are then, respectively, about 2.9 cm 3 s 1 and 2.8 cm 3 s', and are 
essentially the same. The total maximum instantaneous pickup rate is then 5.7 cm 3 s' 1 for the 
nominal source rates assumed here and for this location of Io, which is significantly north of the 
plasma torus centrifugal plane. For the new atomic source rate noted above that are ~4.3 times 
larger, the total instantaneous pickup rate near the exobase would increases to ~25 cm 3 s' 1 , and 
for a finer grid resolution and convergence of the plasma flow lines (not included here), the local 
volumetric rate will increase likely by a factor of several and hence may approach -100 cm' 3 s' 1 . 
This is a significant volumetric source rate at the exobase. It is beginning to approach the 
generic average source rate value of —500 to 1000 cm 3 s 1 that may occur in the Inner Region of 
the denser atmosphere at and below Io's exobase, as simply calculated from a nominal volume 
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integrated Iogenic plasma source of ~1 x 10 28 s' 1 and an atmospheric volume near Io of ~l-2 x 
10 25 era 3 for a spherical atmospheric shell of 200-400 km thickness. 

3.3.2 Ion Pickup Density: Near-Io Scale 

The ion pickup density anywhere in the collection volume about Io can be determined 
from the instantaneous ion pickup production rate. The ion pickup density is obtained by 
integrating the ion pickup production rate along a longitudinal tube (parallel to the streamlines) 
from far upstream to the desired point downstream (see illustration in Figure 4) and dividing this 
integral by the cross-sectional area of the tube and by the local corotational velocity. This 



Figure 4. Plasma Convection Tube for Determination of the Ion Pickup Density. The location of Io 
and several planes perpendicular to Io’s orbital motion are shown in relation to a plasma convection tube 
that extends in the upstream to downstream direction relative to the satellite’s location. 

integration can be expressed alternatively in terms of the time integral of the motion of the ions 
along the tube: 


n 


pickup 


= f P. dt. 

J ton 


tube 
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Assuming curved streamlines that do not depart from the corotational direction, these 
integrations have been performed using the larger volume about lo collected with the coarser 
grid resolution of 2500 km. The ion pickup density is shown in Figure 5 for four selected slices 
in this volume element that are perpendicular to the corotational direction. In Figure 5a, the low 


(30,000 km upstream of lo) 



(7,500 km upstream of lo) 



(at lo) 



(37,500 km downstream of lo) 



RADIAL DISTANCE FROM lO'S ORBIT (R| 0 ) 


Figure 5. Total Iogenic Pickup Number Density for Atomic Oxygen and Sulfur. In a plane 
perpendicular to the undisturbed corotational direction, calculated profiles for the ion pickup 
number density (ions cm 3 ) are shown in (a) 30,000 km upstream lo, in (b) 7,500 km upstream of 
lo, in (c) at lo, and (d) 37500 km downstream of lo. The location of Io's orbit is at the 
intersection of the two dashed lines. Io’s location and size is indicated by the shaded circle. An 
incomplete collisional cascade source at Io’s exobase of 3.2 x 10 27 atoms s ' and 1.6 x 10“ 7 atoms 
s' 1 was assumed, respectively, for oxygen and sulfur. The pickup rates were collected in the 
coarser grid boxes with a spatial resolution of 2500 km. The inner contour in (c) is 20 cm . 
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ion pickup density located 30,000 km upstream of Io can be seen to be asymmetric located about 
Io's orbit with a peak occurring at larger radial distances than Io (positive radial values). This 
pattern occurs primarily because of orbital dynamics that causes the higher-speed neutrals to trail 
Io at radial distances larger than Io's orbit. In Figure 5b, the ion pickup density located 7,500 km 
upstream of Io can be seen to be larger and more symmetric about Io's orbit but still exhibits 
some of the radial asymmetry in Figure 5a because of the contributions of the upstream ions to 
the downstream location. In Figure 5c, the ion pickup density located in a plane containing Io 
can be seen to be significantly larger and more nearly symmetric about the satellite. In Figure 
5d, the ion pickup density located 37,500 km downstream of Io can be seen to be even larger and 
is now fairly symmetric about Io's orbit. This more radially symmetric downstream behavior 
occurs because of the additional source of pickup ions created inside of Io's orbit by lower-speed 
neutrals that, because of orbital dynamics, lead Io and are located radially within Io's orbit. For 
the nominal source rate assumed, the ion pickup density near Io's exobase is larger than 20 cm 
(the inner contour). For a more realistic source that is ~4.3 times larger and for a much smaller 
grid box resolution near Io which should increase the local density by an additional factor of —5- 
10 the ion pickup density near the exobase should then be ~ 400-800 cm' 3 , which is a significant 
fraction of the nominal undisturbed upstream ion density of -2000 cm' 3 . These estimates for the 
ion pickup density near the satellite can be improved by the inclusion of a more accurate 
description for the streamlines near Io and by the inclusion of a description for the Io plasma 
wake in the model calculation of the Iogenic plasma source. 

IV. Summary of Related RCM-J Studies 

New calculations using the Rice Convection Model for Jupiter (RCM-J) incorporating the 
effect of a electromagnetic disturbance at Io to represent a localized plasma source were 
performed and reported at the Fall 2000 AGU meeting (Spiro et al. 2000). These are the first 
calculations that address the question: How does a localized plasma source influence the 
transport convection cells in the Io plasma torus? In these new simulations, the following 
simplifying assumptions were made (although are not required by the model): (1) spin-allowed 
dipole magnetic field, (2) uniform conductivity in Jupiter’s ionosphere, (3) cold torus plasma (kT 
« mQ 2 r 2 /2), and (4) Coriolis force neglected. 
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The moving localized Iogenic plasma source has two major effects: (1) mass loading of 
flux tubes intersecting the source region, and (2) associated pickup currents, which close in the 
Jupiter’s ionosphere, resulting in deflection of the corotational flow around Io. As a first step 
and due to numerical problems, the model does not include the Iogenic plasma source self- 
consistently. The mass is added to a 30° simulation sector every ~13 hours (Io’s period in the 
plasma torus) with a mass addition rate corresponding to 1/12 (307360°) of the normal 10 3 kg/s. 
The pickup current is represented by a 2-D dipole-like current source in the equatorial plane that 
propagates through the 30° simulation sector at ~57 km/s (Io’s orbital speed with respect to the 
corotating frame) once every ~ 1 3 hrs. The size and strength of this dipole current source 
corresponds to complete shielding of the corotational flow from a circular radius of one-fourth of 
a Jupiter radii (~ 10 Io radii), which is physically too large but is a starting point for these initial 
calculations. For simplicity, the radial outflow of the torus is “impounded” by increasing 
Jupiter’s ionospheric conductivity instead of including corotational lag or a ring-current 
impoundment. To increase grid density, only a 60° longitude sector is numerically simulated. 

In the time evolution of the plasma in the 60° longitude sector, tiny ripples are generated 
by Io’s initial passage and because of the interchange instability grow radially outward in finger- 
like regions of enhanced plasma per magnetic flux with a characteristic wavelength that is 
approximately the scale size of the initial electric disturbance. At the second passage of Io, the 
figure-like structures radially beyond Io’s orbit continue to move radially outward, but the base 
of each finger-like structure at Io’s orbital location is undercut and skewed in the direction of 
Io’s motion relative to the outward radial motion of that finger and produces finer and finer 
spatial structures in this region. The plasma inside of Io orbit is relatively undisturbed by the 
passage of Io. Mass added from the second passage of Io moves into the fingers formed in the 
first passage of Io. A similar process occurs on the third passage of Io, with the fingers now 
much longer in their radial extent. 

The inclusion of an Io plasma source makes the physical situation very complicated. The 
simulations, although very interesting at this stage of development, are not sophisticated enough 
yet to produce the details correctly. Some of the general features that may be realistic are: (1) 
the combined effect of the Io-associated pickup current and the interchange-related plasma 
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fingers generate intricate fine structure in the outer part of the main torus near Io’s orbit, (2) 
structure keeps getting finer and finer as Io traverses the modeled region repeatedly, where in the 
real system the fine structure presumably diffuses away when it becomes too fine, (3) the inner 
edge of the torus, which is interchange stable, does not develop the extreme fine structure, and 
(4) each time Io comes around again, new plasma injected into the plasma torus moves outward 
following the path of the pre-existing fingers. Additional calculations (Goldstein et al. 2001) 
including the velocity shear imposed by the corotational lag profiles of the plasma torus near Io’s 
orbit observed by Brown (1994) were presented at the meeting “Jupiter: Planet, Satellites and 
Magnetosphere” held in Boulder, Colorado during June 25-30, 2001. The velocity shear 
significantly alters the development of the radial fingers, producing increasingly complex and 
increasing fine-scale structures down to the smallest scale permitted by the numerical RCM-J 
simulations. Future improvements in the RCM-J include finding a systematic way of diffusing 
the fine structure, including impoundment by corotational lag and/or ring current, and having 
longer simulation times to see if there is an asymptotic behavior for the transport structure. 
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Abstract 


Three-dimensional calculations are presented for the circumplanetary nature of the 
Iogenic plasma source (pickup ions produced by electron and charge exchange processes in the 
plasma torus) created by O and S gases located above Io’s exobase in its corona and escaping 
extended neutral clouds (designated as the “Outer Region”). These calculations are undertaken 
using neutral cloud models for O and S with realistic incomplete collisional cascade source 
velocity distributions at Io's exobase and realistic spacetime loss processes in the plasma torus. 
Two-dimensional distributions produced by integrating the three dimensional information along 
the magnetic field lines are presented for the pickup ion rates, the mass loading rates, the pickup 
conductivity, the pickup radial current, and the net-energy loading rate. On the circumplanetary 
spatial scale, the instantaneous Iogenic plasma source is highly peaked about Io's position on its 
orbit around Jupiter. The degree of orbital asymmetry and its physical origin are discussed and 
overall spatially integrated rates are presented. The relative importance of the Iogenic plasma 
source created by the Outer Region and the Inner Region (region below Io’s exobase) is 
estimated and implications are discussed. 
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1. INTRODUCTION 


Jupiter’s magnetosphere is dominated by heavy-ion plasma (0 + , 0 ++ , S + , S ++ , S +++ ) and has 
its largest density in the Io plasma torus. The plasma torus is located well within the inner 
magnetosphere and is organized about (and extends more radially outward of) Io s orbit, which 
has a radius of 5.9 Jupiter radii (R,). The heavy-ion plasma is supplied by the loss of neutrals 
from Io, which has an atmosphere composed primarily of S0 2 and, to a lesser extent, its detected 
chemically dissociated products (SO, O, S) and trace species (Na, K, Cl). By the Iogenic plasma 
source, we mean the rate of supply of ions, mass, and energy to the plasma torus created by 
neutrals from Io. Although the direct loss of ions from Io’s atmosphere and ionosphere may 
involve photoionization, the Iogenic plasma source is created primarily through the direct and 
more rapid plasma-neutral interactions of electron impact ionization and charge exchange that 
create new pickup ions for the plasma torus. Understanding the three-dimensional nature of the 
Iogenic plasma source is highly relevant to a large number of studies for the Io-Jupiter system 
including the structure and outward transport of the plasma torus and the electrodynamic 
interactions between the plasma torus, Io, and Jupiter. The heavy-ion plasma supplied to the 
plasma torus, for example, forms interesting radial structures, as well as System III longitude and 
dawn-dusk asymmetries, and is ultimately lost by outward radial transport down the tail ot the 
magnetosphere. The creation of the Iogenic plasma source from neutrals from Io can be divided 
into two spatial regions. First, there is a very spatially confined and denser Inner Region 
below Io's exobase where pickup (i.e., both electron impact ionization and charge exchange) ions 
are created and also where Pedersen and Hall conductivities become important in the 
collisionally bound atmosphere and ionosphere. Second, there is a spatially extended and much 
less dense "Outer Region" above the exobase in the bound corona and in the more distant 
escaping neutral clouds where pickup ions are also created. 

The formation of the Iogenic plasma source in the Inner Region is very complex and 
involves the solution of a detailed three-dimensional treatment of the nature of the coupling 
between Io’s atmosphere and its electromagnetic interaction in the magnetosphere. Significant 
efforts to solve this problem have been made in the last few years, but much remains to be 
accomplished. The two-fluid model approach of Saur et al. (1999) explicitly concentrated on 
this Inner Region with the primary purpose of solving the electrodynamic problem for a fixed 
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(unchanging) spherically-symmetric description of the atmosphere. In the one-fluid MHD model 
approach (Linker et al. 1998; Combi et al. 1998; Kabin et al. 2001), the Inner Region and part of 
the near-Io Outer Region are represented in a less detailed fashion by including an empirical 
spatial description for the neutral density in Io’s corona (and more recently in the atmosphere) to 
generate local pickup ions. The Inner Region is also partially represented in the MHD approach 
by the specification of inner boundary conditions, which may be used, for example, to simulate a 
"leaky ionosphere". In contrast to the electrodynamic studies for the Inner Region, two- 
dimensional hydrodynamic studies for the structure, composition and dynamics of Io’s 
atmosphere have been undertaken where the important effects of plasma interactions on shaping 
the atmosphere have been included in simplified manners (Wong and Smyth 2000; Smyth et al. 
2001 ). 


The formation of the Iogenic plasma source in the Outer Region is somewhat less 
complex and can be calculated in three dimensions using neutral cloud models, which have been 
developed over the past two decades (see Section II). These numerical models describe the 
spatial distribution of Io’s neutrals above its exobase and their plasma interactions both in the 
near-Io space of the corona and in the more distant circumplanetary space. Using this approach, 
the general spatial nature of the Iogenic plasma source for the Outer Region was deduced by 
Marconi and Smyth (1996). They argued that the neutral clouds of atoms (or molecules), many 
times mistakenly pictured as being somewhat evenly distributed about Io's circumplanetary orbit, 
should be primarily spatially concentrated in the general vicinity of Io. This must occur 
physically because the typical lifetimes (~1 to 20 hrs) for Io neutrals (H, O, S, SO, S0 2 ) in the 
plasma torus [due to electron impact ionization (or dissociation) and/or charge exchange 
reactions] are short compared to the characteristic time to form a complete toroidal cloud about 
Jupiter at Io's orbit. The destruction of the neutral clouds in the torus should hence produce a 
spatially extended Iogenic plasma source both upstream and downstream of Io that is highly 
peaked at the satellite’s location. 

The Iogenic plasma sources for the Inner Region and Outer Region are coupled in 
interesting ways through various physical processes involving the corotating plasma torus and 
neutrals in the Inner Region. The source velocity distribution of neutrals at Io’s exobase for the 
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Outer Region is created in the Inner Region. This source velocity distribution is composed of 
relatively slow speed neutrals (several km s' 1 , comparable or less than the Io’s surface escape 
speed of 2.56 km s' 1 ) produced primarily by ion-neutral sputtering processes (and quasi- 
collisional thermalization processes) and fast neutrals (10’s to perhaps 100 km s ! ) created by 
charge exchange processes in Io’s local atmosphere. The slow-speed neutrals of the source 
velocity distribution are important dynamically in determining how the neutrals are distributed in 
Io’s corona and in the circumplanetary space near Io’s orbit and consequentially how they 
contribute in spatial distribution and strength to the Iogenic Plasma Source for the Outer Region. 
The high-speed neutrals of the source velocity distribution, although observable for some species 
in the Outer Region, contribute very little to the Iogenic plasma source in the Outer Region but 
are thought to be a major contributor to the Iogenic plasma source in the Inner Region. Hence, 
the study of the rates of energization and escape of neutrals from Io has been used as a direct 
probe to the nature and strength of the local magnetospheric interaction and the Iogenic plasma 
source in the Inner Region (Smyth 1998; Smyth and Marconi 2000). The Iogenic source rates 
from the two regions are thought to be comparable (see Section IV). Understanding the relative 
rates of the Iogenic plasma source in the Inner and Outer Regions is particularly important in 
understanding the structure and energy budget of the plasma torus. The primary density peak in 
the plasma torus, the so-called plasma “ribbon”, for example, is now thought to be formed by the 
Iogenic plasma source created in the Inner Region (Smyth and Marconi 1998a; Hill and Pontius 
1998; Delamere et al. 2001). 

In this paper, the calculations and spatial studies on a circumplanetary scale for the 
Iogenic plasma source created by for atomic oxygen and sulfur in the Outer Region are 
undertaken. In a companion paper (Smyth and Marconi 2001), parallel spatial studies of the 
same Iogenic plasma source on a near-Io scale are presented. The model adopted to calculate the 
Iogenic plasma source in the Outer Region is discussed in Section II. In Section III, the model 
calculations for the circumplanetary distribution of the neutral O and S, the pickup ion rates, the 
mass loading rates, the pickup conductivity, and the pickup radial current are presented. The 
general spacetime nature of the of the Iogenic plasma source is also discussed. In Section IV , the 
model results for the Outer Region are discussed and compared with estimates for the Inner 
Region. A Summary is given in Section V. 
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2. MODEL FOR THE IOGENIC PLASMA SOURCE 


2.1 Neutral Cloud Model 

Three-dimensional calculations for the Iogenic plasma source are performed using the 
neutral cloud model of Smyth and Combi (1988a,b) for satellite corona and extended 
atmospheres. This model, including various refinements, has been applied over the last 25 years 
to studying Io’s Na cloud (Smyth and McElroy 1977, 1978; Smyth 1979, 1983; Pilcher et al. 
1984; Smyth and Combi 1987, 1988b, 1991; Smyth 1992; Smyth and Combi 1997), K cloud 
(Smyth 1992), O cloud (Smyth and Shemansky 1983; Smyth 1992; Scherb and Smyth 1993; 
Smyth and Marconi 2000), S cloud (Smyth 1992), and SO and S0 2 clouds (Scherb and Smyth 
1993; Smyth and Marconi 1998b) as well as the plasma torus and its interactions with Io (Pilcher 
et al. 1984; Smyth and Combi 1987; Smyth 1992, 1998; Smyth and Marconi 1998a, 2000; 
Delamere et al. 2001). The neutral cloud model calculates only the contributions to the Iogenic 
plasma source created by gases above Io's exobase and hence only in the Outer Region. In this 
paper, the Outer Region contributions will be limited to those produced by Io’s O and S neutral 
clouds. In brief, the neutral cloud model calculates a large number of three-dimensional neutral 
particle trajectories by solving the restricted three-body equations of motion, that include the 
satellite and planetary masses, and, if required, solar radiation acceleration. Appropriate 
physically weighted information is computed along the trajectories for each time-step to 
determine the gas density, electron impact ionization and charge exchange rates, and various 
chemical production rates, which may be velocity as well as spacetime dependent (see Smyth 
and Combi 1988a,b for the mathematical formulation and more detailed discussions). By 
following the time history of information along the neutral trajectories, the model calculates the 
instantaneous neutral gas density and consequently the Iogenic plasma source (the instantaneous 
ionization and charge exchange rates) that are produced for gas source conditions specified on 
the exobase of Io moving on its orbit about Jupiter. 

The three-dimensional description of the plasma torus adopted in the neutral cloud model 
is based upon the Voyager 1 plasma properties near western elongation for an offset-tilted dipole 
magnetic field and includes an inherent System III longitude asymmetry and an east-west (i.e., 
local time) asymmetry produced by the presence of an east-west electric field. The inherent 
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System III longitude asymmetry located in the so-called “active sector” (-200-300 System III 
longitude) has a spatial dependence based upon measurements by Morgan (1985a,b) and an 
assumed maximum peak electron density that is a factor of two larger than the normal peak 
plasma torus density outside of the active sector. For the east-west electric field, a nominal value 
is assumed where £=0.025 (the ratio of the east-west electric field to the corotational electric 

field). These Iogenic plasma source rates calculated here hence include the inherently space- 
time dependent lifetime processes created in the corotating Io plasma torus (1) because of its 
oscillatory System III latitudinal motion (~±0.75 Rj) caused by the -7° tilt between the plasma 

torus centrifugal equator plane and Io’s orbital plane and (2) because of its Jovicentric radial 
motions (-0.55 R,) compositely caused in local time by the east west electric field and in System 
III longitude by the effective dipole-offset of the planetary magnetic field. 

2.2 O and S Lifetimes 

The lifetime of oxygen and sulfur atoms in the plasma torus is determined by electron 
impact and charge exchange processes. The reactions and relative importance of these different 
lifetime processes included in the models are summarized in Table 1 for nominal plasma 
properties at Io’s radial distance in the plasma torus centrifugal equator plane at western 
elongation (dusk ansa). The charge exchange cross sections adopted in the neutral oxygen and 
sulfur cloud models are those given by McGrath and Johnson (1989) with refinements provided 
by R.E. Johnson (personal communication, 1989). At Io’s location, charge exchange is seen to 
dominate over electron impact ionization for atomic oxygen while the reverse is true for atomic 
sulfur, which has a lower ionization potential. As noted above, all these processes are inherently 
space-time dependent, and the charge exchange rate is also velocity dependent due to the relative 
motion of the inhomogeneous plasma torus and the atom trajectories. The relative importance 
and absolute lifetime s of processes in Table 1 will thus vary in the plasma torus and depend 
upon local time and the radial, System III latitudinal, and System III longitudinal coordinates. 
As an example, the electron impact (El), composite charge exchange (CE), and total (El + CE) 
lifetimes of atomic oxygen and sulfur are shown, respectively, in Figures 1 and 2 in a 
plasmacentric radial and latitudinal coordinate plane that is located at western elongation (i.e., 
for one local time at the dusk ansa) and at a System III longitude of 230 in the active sector. 
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The motion of Io in the plasmacentric plane at western elongation for all System III longitudes is 
indicated by the gray oval, and Io’s location on the oval for the System III longitude of 230° is 
shown by the dot. For atomic oxygen, the El lifetime (Figure la) has an inner minimum (~40 
hrs) centered on the centrifugal equator near about 5.9 Rj, where the electron density peak 
occurs, and an outer minimum (~40 hrs) centered on the centrifugal equator at ~7 R J? produced 
by the combined effects of a radially outward increasing electron temperature and a decreasing 
electron density. The CE lifetime (Figure lb) for atomic oxygen has only one minimum centered 
on the centrifugal equator at ~5.9 Rj where the ion number densities are largest. The combined 
behavior produces a total O lifetime (Figure lc) at western elongation that has a broader 
minimum (-5.6-7 Rj) with an absolute minimum still determined by CE processes at - 6 R,. For 
atomic sulfur, the El lifetime (Figure 2a) is smaller, as expected, and has a similar spatial 
structure to atomic oxygen with an inner minimums (-7 hrs) and a outer minimum (just below 10 
hrs) centered on the centrifugal equator about -5.9 Rj and 6.8 R J? respectively. The CE lifetime 
(Figure 2b) for atomic sulfur has two relative minimums: a first one (just above 10 hrs) centered 
on the centrifugal equator at about -5.3 R j5 that is created by a local maximum in the S + plasma 
torus density, and a second one (just below 10 hrs) centered on the centrifugal equator about —5.8 
Rj. The total sulfur lifetime (Figure 2c) hence has a deeper minimum centered at -5.9 Rj within 
a broader minimum (-5.6-7 Rj). 

2.3 O and S Source Velocity Distributions at Io’s Exobase 

The source velocity distribution adopted at the exobase for atomic oxygen and sulfur is 
the isotropic incomplete collisional cascade distribution determined for sodium by Smyth and 
Combi (1997). This source distribution peaks at 0.5 km/s, well below the exobase escape speed 
of -2 km/s, and extends in the velocity tail to many tens of km/s. This source distribution 
represents the composite effects of ion sputtering, including heating of the atmosphere. The 
distribution is denoted "incomplete" because the asymptotic tail power-law behavior of E' 4/ \ 
where E is the kinetic energy of a neutral, is overpopulated relative to a complete collisional 
cascade distribution with an asymptotic power-law behavior of E 2 . This overpopulation is 
thought to occur because of the lack of the completion in Io’s atmosphere (where escape occurs 
at the top) of the collision cascade process that moves energy from the velocity tail to the lower- 
velocity core of the distribution. The incomplete collisional cascade flux distribution 
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(j)(v; ot, v b , v M ) expressed in terms of the neutral speed, v, is discussed in detail by Smyth and 
Combi (1988b, in Appendix D) and corresponds to an asymptotic tail power-law behavior given 
by E (a l) and has a peak (most probable) speed of v m that is related nonlinearly to the velocity 
parameters Vb and v M . The core of the source distribution at v=v m provides an effective 
parameter to characterize at Io’s exobase the composite effects of the complex atmospheric 
heating processes that occur below the exobase. For the incomplete collisional cascade process, 
this source distribution should also apply as a first approximation to atomic oxygen and sulfur 
since sodium is a trace species. 

For a set of sodium corona and neutral cloud observations for Io, the incomplete 
collisional cascade distribution was determined to be essentially isotropic at the satellite exobase 
with a peak at v m = 0.5 km/sec (well below the exobase escape speed of -2 km/sec) and a tail 
characterized by a=7/3 that extends to many 10's of km/sec [see Smyth and Combi (1997) for 
the high sensitivity of the fit]. For these parameter values, this normalized source velocity 
distribution is shown in Figure 3 and is to be adopted in this paper for O and S. About 62% of 
the source atoms for this distribution are on non-escape trajectories, while the remaining -32% 
escape Io's gravitational grasp (Smyth 1998). The composite processes below the exobase that 
produce this incomplete cascade distribution are yet to be fully understood from first principles 
by the solution to a complex problem involving the collisional atmosphere and its interaction 
with the plasma torus. The source distribution thus provides a significant source of O and S 
atoms with (1) ballistic trajectories (-62%) that populate only the corona of the satellite, (2) low- 
speed escape trajectories that dynamically are confined relative near Io's orbit, and (3) higher- 
speed escape trajectories that quickly move away from Io’s orbit and contribute primarily to the 
neutral cloud at somewhat larger and smaller radial distances. For sodium the source speed 
distribution is responsible for the creation of the very steep gradient near Io in the spatial 
distribution of sodium within Io's corona as well as the more gradual gradients in the forward 
sodium cloud and the trailing sodium cloud (excluding the higher-speed contributions by the 
directional feature source) near Io’s orbit. Higher- velocity neutral source processes in Io’s local 
atmosphere (i.e., direct ejection and charge exchange) that are responsible for creation of the 
sodium directional feature and the sodium zenocorona are, as would be expected, anisotropic 
with a preferential direction aligned with the plasma corotational motion (see Smyth and Combi 
1997 for overview discussion). For the O and S clouds, these higher-velocity neutral source 


9 



processes will also create neutrals that directly escape Io’s local atmosphere and, along with 
direct-ion escape, occur in the “Inner Region”. These Inner Region contributions to the Iogenic 
plasma source are, however, not the subject of this paper. 

The incomplete collisional source velocity distribution adopted above for O and S is 
assumed to be isotropic at an Io’s exobase radius of 2600 km and have rates determined in the 
following manner. The O source rate for this distribution was determined by Smyth and Marconi 
(2000) using the above Io neutral cloud model from the analysis of two different brightness 
measurements of the O neutral cloud in 6300 A emission. The O source rates obtained in this 
manner were 1.27 x 10 28 and 1.47 x 10 28 atoms s' 1 , respectively, for the measurements of Brown 
(1981) and Thomas (1996). In this paper, an intermediate O source rate of 1.38 x 10 28 atoms s' 1 
(367 kg s' 1 ) is adopted. Since there are currently no direct determinations of the S source rate, a 
value of 0.69 x 10 28 atoms s' 1 (367 kg s' 1 ), one-half of the O atom source rate as would be 
appropriate for a parent S0 2 molecule, is adopted for model calculations in this paper. The total 
source rate at Io’s exobase is then 2.07 x 10 28 atoms s ' (734 kg s ') of which only -38% of the 
atoms in the source distribution are energetic enough to escape Io’s gravity. 

3. MODEL CALCULATIONS 

3.1 Circumplanetary Asymmetric Distributions 

The total nominal lifetime in Table 1 of O and S atoms in the plasma torus centrifugal 
plane at Io’s orbital distance is -18 hrs and 9 hrs, respectively. For O and S atom trajectories in 
the near vicinity of Io’s orbit, the lifetimes are variable and can be several times larger, as 
illustrated in Figures 1 and 2. These O and S lifetimes are, however, short compared to the 
characteristic time of -250 hrs to form a complete 360° toroidal cloud about Jupiter for the slow 
escaping neutrals from Io that populate circumplanetary space near Io's orbit in the source 
distribution of Figure 3. It then follows that the O and S neutral clouds and their Iogenic plasma 
source contributions must be highly peaked at Io's instantaneous position and hence highly 
asymmetric in their spatial distributions about Jupiter. There is a particularly interesting and 
long-standing question to be addressed regarding the nature of this asymmetric circumplanetary 
distribution for the Iogenic plasma source. Near Io, the volumetric rates (ions or mass cm 3 s ') 
will be large but the volume will be relatively small compared to the remaining circumplanetary 
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volume of space. In this large volume of the remaining circumplanetary space, however, the 
volumetric rates will be much smaller. The question then arises as to the overall comparison of 
the total (spatially integrated) rates (ions or mass s ') as a function of circumplanetary distance 
from Io (i.e., is a large volumetric rate times a small volume larger or smaller than a small 
volumetric rate times a large volume?). 

Model calculations for several ion-pickup volumetric quantities have hence been 
undertaken in the proximity of Io's orbit on a complete angular scale and a significant radial scale 
to study the nature of this circumplanetary behavior. For the model description, lifetime 
processes, and source velocity distribution and rates discussed in Section 2, the instantaneous ion 
pickup rates have been calculated in three-dimensions and integrated along the magnetic field 
lines to produce a two-dimensional distribution in the plasma torus with an angular longitudinal 
resolution of 5° and a radial resolution of 0.1 Rj in the radial interval from 4 to 10 Rj. These 
model calculations are performed for Io with a geocentric phase angle of 270° (i.e., western 
elongation or dusk ansa) and an Io System III magnetic longitude of 230° (where Io is north but 
approaching the plasma torus centrifugal equator plane with a magnetic off-set tilted dipole L 
value of 5.88). This location of Io provides a substantial instantaneous Iogenic plasma source 
since the satellite is in the active sector and also the atom trajectories relatively near Io have a 
time history and local-time location such that they interact strongly with the plasma torus. In the 
model calculation presented here, the incomplete collisional cascade distribution is represented 
by nineteen non-uniformly spaced velocity intervals and their corresponding weights. For each 
velocity interval, 1298 atom trajectories are ejected uniformly from Io’s exobase, providing a 
total number of orbits per neutral cloud model calculation of 24,662. The lifetime history for 
each atom trajectory is computed sufficiently far in the past to achieve a steady state description 
for the neutral cloud and the instantaneous Iogenic plasma rates for Io’s location chosen above 
for the model calculation. 

3.2 Distribution of O and S 

It is important to understand the spatial distribution about Jupiter of the O and S densities 
that are shaped by the processes in Table 1 in order to understand the spatial nature of the 
Iogenic plasma source that they create. The O and S neutral densities have been calculated in 
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three-dimensions and integrated perpendicular to Io’s orbit plane to produce a two-dimensional 
distribution with a spatial resolution box size of 2 x 10 4 km [~ 1 1 Io radii (R l0 )]. The model 
calculated O and S neutral column densities viewed perpendicular to Io’s orbital plane are given 
in Figures 4a and 4b, respectively. Both the neutral O and S column density distributions are 
highly peaked about Io’s instantaneous position (small dot at western elongation) and are highly 
asymmetric in longitude. The O and S atoms are concentrated primarily ahead and inside of Io s 
orbit (i.e., forward cloud) where the neutral loss rates in the plasma torus are small because of the 
rapidly radially inward decreasing plasma densities and temperatures. Over a larger 
circumplanetary volume, ranging from near Jupiter to a bit beyond Io’s orbit (and preferentially 
asymmetrically located behind Io’s orbital position), the O and S column densities are a factor of 
~50-100 or so times smaller than in the general vicinity of Io. For both the trailing O and S 
neutral clouds, the column density decreases more rapidly with angular distance from Io than the 
forward cloud because the trailing cloud gases are primarily located outside of Io’s orbit 
(because of orbital dynamics) where the plasma torus loss processes are continuously operative. 
On a circumplanetary scale, the column density for atomic oxygen is significantly larger (factor 
~5) than for atomic sulfur because of the shorter lifetime of atomic sulfur in the plasma torus and 
the factor two smaller sulfur source rate. In Figure 4c, the asymmetry of the column density 
profiles for oxygen and sulfur along the horizontal dashed line are shown in Figures 4a and 4b. 
For O, the peak column density at Io is ~4 x 10 u atoms cm" 2 compared to -5 x 10" atoms cm" 
diametrically opposite Io on the satellite orbit (dashed circle), a factor of about 80 smaller. For 
S, the peak column density at Io is ~1.8 x 10 13 atoms cm 2 compared to ~8 x IO 10 atoms cm' 
diametrically opposite Io on the satellite orbit (dashed circle), a factor of about 225 smaller, 
reflecting the smaller lifetime for S. 

3.3 Distribution of the Pickup Ion Rates 

The loss processes for the neutrals in Table 1 are the same processes that create the 
pickup ions for the plasma source and hence determine the Iogenic plasma source for the Outer 
Region. In discussing pickup ion rates, it is important to differentiate the “new pickup ion rate” 
from the “net pickup ion rate”. All of the reactions in Table 1 for electron impact (El) ionization 
and charge exchange (CE) processes create new ions and hence contribute to the “new pickup 
ion rate”. These newly created ions, in the presence of the magnetospheric local electric and 
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magnetic fields, create pickup currents as they are accelerated to the local plasma corotational 
speed and are energize gyrationally. All of the El reactions, but only some of the CE reactions in 
Table 1, however, create net ions for the plasma torus and hence are able to contribute to the “net 
pickup ion rate”. Symmetric charge exchange reactions (e.g. S + S + — > S + + S ; O + 0 ++ — » 0 ++ 
+ O) do not contribute to the net pickup ion rate or to the corresponding net-mass pickup ion- 
loading rate of the plasma torus. Asymmetric charge exchange reactions (e.g. S + 0 + — > S + + O ; 
O + S + — > 0 + + S) do not contribute to the net pickup ion rate but can contribute either 
positively or negatively to the net-mass pickup ion-loading rate. Finally, charge transfer 
reactions (e.g., O + S ++ — » 0 + + S + ) contribute positively both to the net pickup ion rate and the 
net-mass pickup ion loading rate. Since the symmetric charge exchange reactions are the 
dominant CE processes in Table 1 , it follows that the CE processes play a much less important 
role compared to the El processes in their contributions to net pickup ion rate and to the net-mass 
pickup ion loading rate. The CE processes, however, provide a strong contribution to the pickup 
ion conductivity and current and also to the net-energy loading rate for the plasma torus. 

The instantaneous new ion pickup rates produced by electron impact (El) ionization and 
the set of charge exchange (CE) processes in Table 1 are shown, respectively, for the O neutrals 
in the Outer Region in Figures 5a and 5b and for the S neutrals in Figures 5d and 5e. One- 
dimensional profiles for the new pickup ion rates for the El and CE processes are also shown 
along an east-west dashed line (i.e., the line through Io and Jupiter) in Figures 5c for O and in 
Figure 5f for S to better display their circumplanetary asymmetry. For both oxygen and sulfur, 
the new ion pickup rates are highly peaked about Io (small dot at western elongation). An 
examination of the nature of the pickup ion distributions near Io at a higher-spatial resolution 
than considered here is presented in a companion paper (Part II). All of the pickup ion rates in 
Figure 5 decrease rapidly inside of Io’s orbit because the plasma torus densities and temperatures 
decrease rapidly there with decreasing radius. 

A comparison of the new pickup ion rates for El created by O (in Figure 5a) and S (in 
Figure 5d) shows that they have almost identical circumplanetary distributions near and inside of 
Io's orbit. This same morphology occurs because the factor of ~5 larger column density for O 
than for S (in Figure 4) is compensated by the factor of ~5 times smaller electron ionization 
lifetime for S than for O. Well outside of Io's orbit, the new pickup ion rate produced by O 
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becomes larger than that produced by S, since the factor of ~5 enhancement in the O column 
density is no longer fully compensated by the electron ionization lifetime that, at these larger 
distances (see Figures 3a and 4a), is now only a factor of ~2 times smaller for S than O. 
Examination of the one-dimensional profiles peaks in Figures 5c and 5f, which occur at Io and 
diametrically opposite Io on its orbit, shows that the new pickup ion rate for El pickup processes 
in Figures 5a and 5d differ at the two positions by a factor of -600 for O and -1000 for S. 

For CE processes, the pickup ion rates created by O (Figure 5b) and for S (Figure 5e) are 
highly peaked in a narrow radial interval occupying angularly incomplete forward and trailing 
arms that reach out from Io and are located from just outside to just inside of Io's orbit, where the 
charge exchange lifetimes have their minimum values. In these arms, the new pickup ion rate 
created by O is larger than for S by a factor -4, which is caused by the factor of -5 larger O 
column density and a slightly smaller S electron charge exchange. At radial distances well 
outside Io’s orbit, the new pickup ion rate for O relative to S becomes even larger (-20) because 
the S column density decreases more rapidly with Jovicentric distance than the O column 
density. Examination of Figures 5c and 5f shows that the ratio of the new pickup ion rate at Io to 
that diametrically opposite Io on its orbit is -340 for S and -140 for O. 


The spatially integrated new pickup ion rates produced by the El and CE processes in 
Figure 5 and also the associated net pickup ion rates are summarized in Table 2. For El 
processes, the integrated new pickup rates for O (1.28 x IO 27 ions s' 1 ) and S (1.75 x 10 27 ions s ‘) 
are comparable, as expected, with the larger rate for the S with the shorter ionization lifetime and 
hence the larger source rate near Io (see Figures 5c and 5f). For CE processes, the integrated 
new pickup rate for O (3.47 x 10 27 ions s' 1 ) is considerably larger than for S (1.05 x 10 27 ions s' 1 ) 
and reflects the shorter charge exchange lifetime and larger source for O near Io. For O and S 
combined, the El (3.03 x IO 27 ions s' 1 ) and the CE (4.52 x 10 27 ions s' 1 ) contributions are 
comparable and are -40 and 60%, respectively, of the overall integrated new pickup ion rate of 
7.55 x 10 27 ions s' 1 . The spatially-integrated net pickup ion rate (3.85 x IO 27 ions s' 1 ) is about 
one-half this rate and is created mainly (79%) by an El contribution (3.03 x 10 27 ions s ') with a 
smaller (21%) CE contribution (0.82 x 10 27 ions s' 1 ). 
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3.4 Distribution of the Total-Mass and Net-Mass Loading Rates 

The instantaneous equivalent (El + CE) total-mass loading rate (M (EI) + M <CE> ) for the Io 
plasma torus is determined directly from the pickup ion rates for El and CE processes, denoted 
here by S EI and S EE for O and by Sf and Sf E for S, respectively: 

M <e,) = m 0 S EI + m s S s EI da) 

M (CE) = m c Sf + m s Sf E (lb) 

where m 0 and m s are the mass of an O and S ion, respectively. The instantaneous net-mass 

loading rates M (ne,) are determined straightforwardly by adding the rate of creation of ion mass 
for all El processes and of the net ion mass contributions for all CE processes (see Smyth and 
Combi 1988a for a more detailed formulation). The circumplanetary distribution of the 
instantaneous net-mass loading rates and the equivalent total-mass loading rate integrated along 
the magnetic field line are shown, respectively, in Figures 6 and 7. 

The nature of the angular distribution of the net-mass loading rates about Jupiter has been 
an interesting and open question. Suggestions for the shape of the distribution have ranged from 
the extremes of an angular distribution that is uniform in a toroidal volume about Jupiter to a 
distribution that is very highly localized in the immediate vicinity of Io's instantaneous position. 
The distribution in Figure 6a has the character of an immediate state: the source is both very 
highly peaked at Io's instantaneous position and also exists at much lower levels all around 
Jupiter. To quantify this angular variation more succinctly, the two-dimensional distribution in 
Figure 6a is integrated radially within 5° pie-shaped angular intervals to produce the average 
angular distribution shown in Figure 6b. Examination of this angular distribution shows that it 
may be effectively divided into a highly peaked component near Io that is located approximately 
within ±20° of Io's instantaneous position and a more gradually decreasing component that has a 
minimum located approximately on the opposite side of Jupiter from Io. The steeper slope, 
occurring for System III longitudes that are larger than Io, corresponds to the plasma source that 
is located behind (i.e., upstream) of Io and primarily outside its orbit in the hotter region of the 
plasma torus. The less steep slope for System III longitudes smaller than Io corresponds to the 
plasma source that is ahead (downstream) of Io and primarily within (but near) its orbit, where 
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the neutral loss processes are still rapid enough to contribute significantly to the Iogenic source. 
The total (El + CE) pickup ion equivalent mass-loading rate has a similar morphology and is 
shown in Figures 7a and 7b. 

The spatially integrated net-mass loading and total equivalent mass-loading rates for 
Figures 6 and 7 are, respectively, 154 and 277 kg s' 1 . In Table 2, these “Outer Region ’ 
integrated rates are separated into their contributions from electron impact ionization and charge 
exchange. The net-mass loading rate of 154 kg s’ 1 is created primarily by electron impact 
ionization (127 kg s' 1 ) but also has a small contribution from charge exchange processes (27 kg s' 
') because of asymmetric charge exchange and charge transfer reactions in Table 1. The net- 
mass loading rate within the ±20° of the Io peak in Figure 6 is about 109 kg s' 1 (71% of the 
spatially integrated rate) with only 45 kg/s (29%) located within the remaining 320° angular 
interval about Jupiter. The total mass-loading rate within the ±20° of the Io peak in Figure 7 is 
about 197 kg/s (also 71% of the spatially integrated rate) with only 80 kg/s (29%) located within 
the remaining 320° angular interval about Jupiter. 

3.5 Distribution of the Net-Mass per unit Magnetic Flux Source Rate 

One of the volumetric quantities of interest for transport studies in the plasma torus, such 
as those undertaken with the Rice Convection Model for Jupiter (RCM-J; see, for example, 
Pontius et al. 1998), is the source rate for the net-mass per unit magnetic flux integrated along 
the magnetic field, which is denoted by T| : 

f M (ne,) , 

ri = ds (2) 

1 J B 

where B is the local value of the magnetic field, and “ds” is the differential distance element 
along the magnetic field line. For the combined O and S pickup ion rates discussed above, the 
two-dimensional spatial distribution for this instantaneous source rate T| is shown in Figure 8. 

The circumplanetary distribution for the net-mass per unit magnetic flux in Figure 8a is 
highly peaked at Io’s location and is enhanced in both forward and trailing arms near Io’s orbit 
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and in an extended radial region outside lo’s orbit and trailing Io’s orbital location. Upon 
comparison with the net-mass distribution in Figure 6a, the distribution in Figure 8a is seen to be 
skewed by the ~V nature of the factor of 1/B in the definition for T| (where B= 1835 nT near 
Io’s orbit) so as to be relatively lower by about a factor of < 2 inside of Io’s orbit and to be 
relatively enhanced by up to a factor of 4 or more outside of Io’s orbit. In Figure 8b, the degree 
of asymmetry of the distribution is shown in a profile taken along the dashed horizontal line in 
Figure 8a through Io. The maximum value at Io is ~6 x 10' 8 kg m 2 s' 1 T 1 and is a factor of ~ 600 

times larger than the peak diametrically opposite Io on its orbit. 


3.6 Distribution of the Pickup Conductivity, Current, and Power 

The pickup ion rates produced by the O and S neutrals also produce a pickup 
conductivity G pu through the introduction of ion-electron pairs, which are separated on average 

by their gyration radii at the initial times of their creation (Goertz 1980). The pickup 
conductivity depends upon all El and CE processes and is defined by 


M ,EI) +M (CE) 

G pu ~ B 2 

The height-integrated pickup conductivity Z pu is defined by integrating 
field line 


(3) 

G along the magnetic 


I 


pu 



(4) 


In an Io frame with the corotational radial electric field denoted by E , the pickup conductivity 
generates a radial pickup current density J pu (amps m‘ 2 ), which is given by 


J 


pu 



E . 


(5) 


The instantaneous radial pickup current I pu (amps) at a function of radius and System III 
longitude defined through a surface with a differential area dA = £ds, where one dimension is 
along the field-aligned height-integrated direction “s” and the second dimension is along the 
System III longitude direction with a arc length f corresponding to the adopted circumplanetary 
grid resolution of 5°, can also be determined from J pu 
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The power density P pu (watts m' 2 ) required to produce the new pickup ion rates in Io’s reference 
frame is determined from J pu as follows: 


P = f J -E ds = fa 

pu J P u J pu 


E 2 ds 


(7a) 


and may be written alternatively as 

P pu = i j(M (E,) + M (CE) )V c 2 or ds + - J(M (E,) + M (CE> )V c 2 or ds . (7b) 

2 2 

where the local corotational pickup speed V cor is 



ExB _ E 
B 2 ” B ' 


( 8 ) 


The expression (7b) for the power density P pu , which has been divided into two equal parts, has 

an obvious interpretation. The first half of the power density is the work required by the electric 
field to accelerate the center-of-mass motion of the pickup ions to the local corotational speed of 
the plasma torus, while the second half is the power required to energize the gyrational motion of 
the ions relative to their center-of-mass motion. This center-of-mass motion of the pickup ions is 
reflected in the expression for the force exerted on the pickup ions 


J u xB = G E x B = (M |E,) + M ,CE) )^ = (M ,EI) + M (CE) )V cor . (9) 

pu P u B 

The gyration power of the new pickup ion energy is the rate of energy at which new ions are 
created to heat the plasma torus, i.e., the “new pickup ion energy loading rate”. It is, however, 
larger than the “net pickup ion energy loading rate” for the plasma torus since symmetric and 
asymmetric charge exchange reactions transform old cooler ions of the plasma torus into fast 
neutrals that rapidly escape from the plasma torus with the energy of the old cooler ions. 

The circumplanetary distribution of the height-integrated pickup conductivity £ pu is 
shown in Figure 9a. The pickup conductivity Z pu is seen to be highly peaked at Io’s location 


and is enhanced in both forward and trailing arms near Io’s orbit and in an extended radial region 
outside Io’s orbit and trailing Io’s orbital location. Upon comparison with the total-mass 
distribution in Figure 7a, the distribution in Figure 9a is seen to be heavily skewed by the ~L 6 
nature of the factor of B 2 in the definition for G pu so as to decrease more rapidly inside of Io’s 

orbit and to decrease much more slowly outside of Io’s orbit. In Figure 9b, the degree of 
asymmetry of the distribution is shown in a profile taken along the dashed horizontal line in 
Figure 9a. The maximum value at Io is ~2.5 x 10' 3 mho and is a factor of ~ 450 times larger than 

the peak diametrically opposite Io on its orbit. 

The circumplanetary distribution of the radial pickup current I pu is shown in Figure 10a. 

The radial pickup is highly peaked at Io’s location and asymmetrically distributed about Jupiter 
with enhanced forward and trailing arms that extend further outside of Io’s orbit near the satellite 
and become confined very near Io’s orbit at much larger angular distances from the satellite. In 
Figure 10b, the peak current at Io is ~4 x IO 4 amps while that at the location opposite Io on its 
orbit is -100 amps, an asymmetric factor of 400 smaller. 

The circumplanetary distribution of the net-energy loading rate is shown in Figure 11a 
and is highly peaked at Io and has a morphology in asymmetric angular similar to the radial 
pickup current in Figure 10a. In Figure lib, the peak net-energy loading rate at Io is - 6 x 10 4 
watts m 2 as compared to -1.8 x 10' 6 watts m' 2 diametrically opposite Io on the satellite orbit, an 
asymmetry factor of over 300. The spatially integrated net-energy loading rate for Figure 1 la is 
1.03 x 10' 2 watts and is dominated by the contribution from the large peak at Io. 

3.7 Spacetime Variability of the Iogenic Plasma Source 

To understand the general nature of the space-time variability of the Iogenic plasma 
source rates (e.g., the net ion-loading rate, net mass-loading rate, equivalent mass loading-rate 
for momentum loading, and energy-loading rate), it is important to understand the relationship 
between the instantaneous neutral source rate at Io, the neutral cloud instantaneous accumulation 
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rate in the circumplanetary space, and the instantaneous loss rates of neutrals in the 
magnetosphere, the latter of which produces the Iogenic plasma source. The relationship among 
these three quantities has been derived mathematically using the general "ensemble packet" 
formalism developed by Smyth and Combi (1988a). The results are presented briefly below. 

Consider on Io's exobase a solid angle element d£2 at the location Q and atoms ejected 
from this element with a velocity in the range to to to + dco . We denote the neutral source rate 

A ^ 

(atoms per second) originating with the conditions (to, Q) at time T by <{>((0, £2, T) dco dQ where 
0(0), f2,T) is the flux per unit d£2 and per unit dco. We denote by N(to, Q, T) dtb dQ the 
number of neutrals at time T along the complete neutral trajectory formed prior to the time T and 
originating with the conditions (cb/Q). We furthermore denote the neutral cloud instantaneous 
accumulation rate at time T for atoms originating with the conditions (to, Q) for past times up to 
time T by 


dN(cb / fl,T) 
” dT 


dcodQ. 


( 10 ) 


Finally, we denote by L(co, £2, T) dco d£2 the instantaneous neutral cloud loss rate at time T 
caused by magnetospheric processes for the complete trajectory originating with the conditions 
(to, Q). The desired relationship between these three quantities may then be expressed as 
follows 


m o, T) - ♦<““>(», n, T) = dN<0 l n '9 + L(S, 6 , T) (11) 

dT 

where <j/ ex,t *(co, £2,T) is the instantaneous atom flux along the trajectory that is lost from the 
trajectory because of collision with Io, collision with Jupiter, or that remains undepleted and 
unchanged on the trajectory having been ejected from Io at times prior to time T. 

The equation (11) is a conservation statement at time T for the complete trajectory 
originating with the conditions (63,0). Since it is valid for all neutral trajectories from Io s 


20 



exobase, it also follows that it is true for the neutral cloud as a spatially integrated whole by 
simply integrating over all conditions at the exobase. The left hand side of equation (11) is the 
instantaneous net flux at time T from Io that is available to either increase the neutral cloud 
abundance along the entire trajectory or to supply the neutral cloud loss rate caused by the 
magnetospheric processes. Note that this neutral cloud loss rate may be smaller than the net flux 
if the abundance of neutrals along the trajectory is increasing or may be larger than the net flux if 
the abundance of neutrals along the trajectory is decreasing. Hence, the neutral cloud is a 
reservoir for the neutrals that can increase in abundance when the global loss rate becomes 
smaller and can decrease in abundance when the global loss rate becomes larger. The time 
variability in the loss rate is caused by the oscillation of the plasma torus about Io’s orbital plane, 
the System III longitude asymmetries of the plasma torus, and the east-west asymmetries of the 
plasma torus. Hence, the instantaneous strength of the Iogenic plasma source determined by the 
neutral cloud loss rate may be larger than the instantaneous net flux at time T from Io if the 
neutral cloud loss rate is decreasing and conversely may be smaller than the instantaneous net 
flux at time T from Io if the neutral cloud loss rate is increasing. Hence, it is not the 
instantaneous net flux at time T from Io that determines the instantaneous strength of the Iogenic 
plasma source but the complete time history at time T of the neutral cloud. 

4. DISCUSSION 

4. 1 Comparison of Inner and Outer Region Iogenic Sources 

Several estimated rates for the spatially integrated Iogenic plasma source are summarized 
in Table 4. These rates combine in various ways contributions from the Inner Region and the 
Outer Region and hence can be considered to be approximate estimates for the total Iogenic 
plasma source. The estimates of Smyth (1992) and Smyth and Marconi (2000) were obtained 
from neutral exobase source rates for Io determined by modeling the spatial brightness 
distribution of the neutral clouds and from additional scaling arguments to extend the velocity 
range to include high-speed charge exchange processes. The estimates of Smyth (1992) used 
monoenergetic source conditions to describe the slower neutrals at Io’s exobase rather than the 
more accurate description of the incomplete collisional cascade distributions used here and in 
Smyth and Marconi (2000). Hence, the more recent estimates must be considered to be more 
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accurate. The estimates of Smyth and Marconi (2000) include rates for the Outer Region, 
essentially those presented in this paper, and additional estimated rates for the Inner region that 
are included only for charge exchange processes. The estimates of Bagenal (1997) were 
obtained from examination of the Galileo PWS electron data (Gurnett et al. 1996) and PLS 
plasma data (Frank et al. 1996). The rates for Combi et al. (1998), Linker et al. (1998), and 
Kabin et al. (2001) were obtained from MHD modeling of the Galileo PLS plasma (Frank et al. 
1996) and MAG magnetic field data (Kivelson et al. 1996) and are listed for different scenarios, 
which provided marginal fits for most of the data but failed to reproduce the two local minima in 
the magnetic field on the wake flanks. Combi et al. (1998), who did not solve for the magnetic 
field inside of the satellite, considered two different inner boundary conditions [fixed (f) and 
reflective (r), in Table 4] at an altitude of 150 km to span extreme conditions for the particle and 
magnetic field interactions with Io and its ionosphere. For the "reflective" boundary conditions, 
Io and its ionosphere are treated as a perfect conducting impenetrable sphere so that both the 
magnetic field and plasma flow are reflected. For the "fixed" boundary conditions, the magnetic 
field and particles can penetrate the ionospheric boundary and are treated at the boundary by 
setting the fluid velocity to zero, the density to ionospheric values, and fixing the magnetic field 
so that only a small expected change occurs in its value inside Io’s solid body. Linker et al. 
(1998), who solved for the magnetic field inside of Io, considered the case of an Iogenic plasma 
source with both a magnetized and non-magnetized Io [labeled, respectively, (m) and (n) in 
Table 4], Kabin et al. (2001) added to the earlier treatment of Combi et al. (1998) new pickup 
ion terms into the equations and included, in addition to a power law description for the neutral 
density in Io’s corona, an exponential scale-height term to describe Io’s atmosphere nearer the 
surface and also a day-night asymmetry in the neutral density distribution. 

The estimates of Smyth (1992) and Smyth and Marconi (2000) have a ratio of the charge 
exchange rate to the net mass-loading rate of ~1.6 and -3.3, respectively, which is similar to the 
values of -1.1, -1.8, and 3, respectively, for Kabin et al. (2001), Combi et al. (1998), and Linker 
et al. (1998) but is in striking contrast to the values of Bagenal (1997), who concluded that there 
was no significant charge exchange rate near Io at the time of the Galileo flyby. Since charge 
exchange reactions are always operative and physically dominant over electron ionization rates 
at Io's position in the undisturbed plasma torus, this conclusion is suspect. The conclusion of 
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Bagenal (1997) can be traced effectively to an arbitrary division of the flux of the "new ionized 
ions" from the "background ions" in the PLS and PWS measurements. In Table 4, the estimates 
for the total new pickup ion rate of 5 x 10 28 ions s' 1 for the more extreme reflective boundary 
condition model [labeled (r)] of Combi et al. (1998) is likely too large. This larger source rate 
was required to supply the same amount of plasma density to the plasma wake region, where in 
this case both the particles and the magnetic field are highly reflected by the inner boundary and 
are therefore distributed over a larger spatial volume. The smaller total source rate (n) of Linker 
et al. (1998) is also likely too small since it was determined for the case of an assumed intrinsic 
magnetic field for Io, which does not now appear to be supported by MHD modeling of more 
recent Galileo flyby magnetospheric data (Linker et al. 2001). With these limitations then, the 
best estimates obtained from Table 4, which represent the combined contributions of the Inner 
and Outer Regions, are very similar and are in the range -1.7-2.7 x 10" 8 ions s' 1 for the total new 
pickup ion rate and -600-980 kg s' 1 for the total mass loading rate. 

The estimates of the new mass loading rate for the Iogenic plasma source for the Inner 
Region determined in this paper and for the combined Inner and Outer Regions obtained from 
Table 4 are summarized in Table 5. By subtracting these rates, a rough estimate for the new 
mass loading rate for the Inner Region is therefore -300 to 700 kg s' 1 so that the rate for the Inner 
Region is -1 to 2.5 larger that for the Outer Region. The Inner Region is, however, highly 
concentrated in a small volume around Io. In contrast, the primary peak of the Outer Region at 
Io, which contains -71% (125 kg s' 1 ) of the new mass loading rate, is located within a large 
volume with angular dimension of ±20° longitude (i.e., ±81 R IO ) about Jupiter. A more detailed 
description of the spatial distribution of Iogenic plasma source for the Outer Region within this 
large volume is the subject of the companion paper. Hence volumetrically, the Iogenic plasma 
source near Io must be dominated by the Inner Region contribution. Similarly, the net energy 
pickup rate of -1 x 10 12 watts calculated for the Outer Region is less than the -2-5 x 10 12 watts 
estimated for the total energy budget of the plasma torus. The Inner Region should thus also be 
expected to have a comparable or larger contribution, which should balance the energy budget. 
The asymmetries in the Iogenic plasma sources at Io will therefore be much larger than those 
presented in Section 3 for the Outer Region source only. Given this perspective, it should not be 
so surprising that the plasma torus ribbon, the dominant radial density feature in the plasma 
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torus, is created by the Iogenic plasma source created in the Inner region (Smyth and Marconi 
1998a; Hill and Pontius 1998; Delamere et al. 2001) and that the energy budget of the plasma 
torus is very dependent upon the power ribbon created by this Inner Region source (Smyth and 
Marconi 2000). 

4.2 Additional Iogenic Sources 

The Iogenic plasma source calculations presented in Section 3 were only for pickup ions 
created by O and S gases in the Outer Region. These rates may only be a lower limit if other 
species (S0 2 , SO, O,, S 2 ) in Io’s atmosphere can also contribute significantly to the Iogenic 
plasma source in the Outer Region. The abundances of 0 2 formed through atmospheric 
chemistry and S, released by volcanoes are both thought to be too small on a global scale to 
provide a significant source of pickup ions. The most abundant species S0 2 and the next most 
abundant species SO (-10% of S0 2 ) in Io’s atmosphere are expected to undergo rapid electron 
impact dissociation in Io’s upper atmosphere (Smyth et al. 2001) producing primarily O and S 
exothermic atomic products. The electron impact production rate for molecular ions is small 
compared to that for molecular neutral products. Much of the exothermic O and S source rates at 
the exobase may, however, already effectively be included in the source rates for the Outer 
Region adopted here, which are calibrated by observations for atomic oxygen. The direct escape 
of S0 2 and SO from the exobase may, however, provide some spatially extended source of 
exothermic O and S to the circumplanetary space, although SO is likely to be less important 
because of its order of magnitude smaller abundance relative to S0 2 in Io’s atmosphere. The SO, 
and SO exobase source rates are not well known, and their exothermic O and S energy 
distributions have not been measured. An S0 2 exobase source rate of 4 x 10 27 molecules s 1 was 
given by Smyth and Marconi (1998b) based upon the analysis of ion cyclotron waves for a 
maximum SO, + density profile determined by Huddleston et al. (1997) and the same incomplete 
collisional source distribution adopted in this paper for O and S. Hence only -38 percent of the 
SO, molecules are on escape orbits, providing an SO, escape flux of -1.5 x 10 27 molecules s' 1 . 
From a second analysis of the same SO, + density profile by Crary and Bagenal (2000) using a 
simpler monoenergetic velocity of 2.48 km s’ 1 for the source of SO, at Io’s exobase, a similar 
SO, escape rate of a 1-3.5 x 10 27 molecules s' 1 was suggested. If the S0 2 exobase source rate is 
indeed this large and if the exothermic products are not very energetic so as not to dilute their 
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local circumplanetary density contributions, then S0 2 may possibly provide an additional 
significant contribution to the Iogenic plasma source in the Outer Region. Information regarding 
the electron impact cross-sections and their exothermic energy distributions for the products in 
the molecular dissociation of S0 2 (and SO) is needed to address this question. 

5. SUMMARY 

The instantaneous pickup ion rates for the Iogenic plasma source created O and S in the 
Outer Region through electron impact and charge exchange processes in the plasma torus have 
been calculated in three dimensions on a circumplanetary scale for one location of Io in System 
III longitude and local time. The pickup ion rates for O and S are presented in two-dimensional 
distributions by integrating the three-dimensional rates along the magnetic field line. Similar 
two-dimensional distributions for a number of volumetric quantities determined from these 
pickup ion rates are presented for the net-mass and total mass loading rates, the mass per unit 
magnetic flux rate, the pickup conductivity, the radial pickup current, and the net-energy loading 
rate for the plasma torus. All of the two-distributional distributions are highly peaked at Io’s 
location and hence highly asymmetric about Jupiter. The physical origin and the amount of 
asymmetry for the volumetric quantities are discussed. The spatially integrated rates for the net- 
mass and total mass loading rates are 154 kg s' 1 and 277 kg s' 1 , respectively, with 71% of the rate 
confined to a very large peak located within ±20° longitude of Io’s orbital location. The spatially 
integrated rate for the mass per unit magnetic flux is 1.1 x 10 8 kg s 1 T ‘. Maximum values for 
the pickup conductivity and the radial pickup current at Io’s location are 2.5 x 10 3 mho and 4 x 
10 4 amps. The spatially integrated net-energy loading rate is 1.03 x 10 12 watts and is dominated 
by the contributions from the peak centered on Io. 

The space-time variability of the Iogenic plasma source is discussed and is shown to 
depend upon the complete time history of the neutral clouds from which they are created. The 
Iogenic plasma source calculated here for the Outer Region is compared to various estimates for 
the Iogenic plasma source created in the combined Inner and Outer Regions, and the mass source 
rate for the Inner Region (-300-700 kg s ') is estimated to be -1 to 2.5 times larger than that 
calculated for the Outer Region (277 kg s '). Other possible sources of pickup ions for the Outer 
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Region are also discussed. The Inner Region is therefore expected to be the dominant 
contributor of ions to the formation of the plasma torus ribbon, the maximum density peak in the 
radial structure of the plasma torus, and a comparable if not larger contributor to the energy 
budget of the plasma torus compared to the Outer Region. Future improvements in calculating 
the Iogenic plasma source include the addition of near-Io plasma effects in the satellite flank and 
plasma wake and the more complete specification of source rates for neutrals at Io’s exobase. 
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Table 1 . O and S Loss Processes in the Plasma Torus near Io’s Orbit at Western-Elongation 





Relative Rate, 
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40000 
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38000 
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100.00 


Atomic Sulfur 
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Table 2. Volume Integrated Pickup Ion Rates 


Neutral Species 

Electron Impact Ionization 
(10 27 ions s' 1 ) 

Charge Exchange 
(10 27 ions s' 1 ) 

Total 

(10 27 ions s' 1 ) 

New Pickup Ion Rates 
O 

1.28 

3.47 

4.75 

S 

1.75 

1.05 

2.80 

O + S 

3.03 

4.52 

7.55 

Net Pickup Ion Rates 

3.03 

0.82 

3.85 
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Table 3. Net and Total Mass Loading Rates’* 


Magnetospheric Process Net Mass Loading Rate Total Mass Loading Rate 

(kg/s) (_%T (kg/s) (%) a 


Electron Impact Ionization 

127 

17.3 

127 

17.3 

Charge Exchange 

27 

3.7 

150 

20.4 

Total: 

154 

21.0 

277 

37.7 


a Based upon the combined O and S exobase upward source rate of 734 kg/s (100%); about 62% of the source rate 
produces ballistic atom trajectories that are confined to Io's corona. 
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Table 4. Estimates for the Total Iogenic Plasma Source 


Net Mass-Loading Charge Exchange Total 


(10 28 s ’) 

(kg s' 1 ) 

(IO 28 s 

') (kgs' 1 ) 

'c/2 

oc 

r-i 

o 

(kg s ') 

Reference 

Neutral Clouds 






0.57 

225 

1.15 

367 

1.72 

592 

Smyth (1992) 


154 


510 

1.95 

664 

Smyth and Marconi (2000) 

Galileo Flvbv 







0.53-1.7 

188-600 

0 

0 

0.53-1.7 

188-600 

Bagenal (1997) 

MHD Simulations 






0.8 

300 (f) 

1.5 

540 (f) 

2.3 

840(f) 

Combi et al. (1998)“ 

1.8 

650 (r) 

3.2 

1150 (r) 

5.0 

1500 (r) 


0.5 

200 (n) 

1.5 

600 (n) 

2.0 

800 (n) 

Linker et al. (1998) b 

0.3 

120 (m) 

0.9 

360 (m) 

1.2 

480 (m) 


1.3 

470 

1.4 

510 

2.7 

980 

Kabin et al. (2001) 


a The two cases are for fixed (0 and reflective (r) inner boundary conditions about Io. 
b The two cases are for no (n) intrinsic magnetic field and an assumed (m) intrinsic magnetic field for Io. 
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Table 5. Spatial Distribution: Total Ion Pickup Loading Rate 

Mass Rate 

(kg/s) 

Iogenic Plasma Source (Outer Region) 

Exobase to the ±20° Base of the Io Peak 197 

Remaining 320° of System III Longitude about Jupiter 80 

Total Outer Region Source 277 

Iogenic Plasma Source (Inner and Outer Regions) 600 -980 
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FIGURE CAPTIONS 


Figure 1. Atomic Oxygen Lifetime in the Plasma Torus at Western Elongation. The 

electron impact (El), composite charge exchange (CE), and total (El + CE) lifetimes, in units of 
hours, of atomic oxygen are shown in (a), (b), and (c), respectively, in a plasma torus plane 
located at a local time corresponding to western elongation (dusk ansa) and at a System III 
longitude of 230°. The radial ( L) and vertical (z) plasmacentric coordinates (Smyth and Combi 
1988b) are defined so that the density structure of the plasma are fixed in location and Io moves 
around the gray oval as a function of System III longitude and at 230° is located at the position of 
the gray dot. 

Figure 2. Atomic Sulfur Lifetime in the Plasma Torus at Western Elongation. Same caption 
as Fig. Caption 2, except for atomic sulfur. 

Figure 3. Incomplete Collisional Cascade Source Velocity Distribution. The source speed 
distribution (normalized to unit area) adopted for the flux of O and S at Io’s exobase is that 
determined for Io’s sodium corona and neutral clouds by Smyth and Combi (1997). The 
distribution peaks at 0.5 km s' 1 and has an incomplete collisional cascade tail characterized by 
the parameter a = 7/3 (see text ; also Smyth and Combi 1988b). 

Figure 4. Distribution of Atomic Oxygen and Sulfur about Jupiter. For a line of sight 
perpendicular to Io’s orbital plane, the model calculated circumplanetary distribution of the 
column density (atoms cm' 2 ) is shown for atomic oxygen in (a) and for atomic sulfur in (b). The 
relative location of Jupiter (half lit disk), Io’s orbit (dashed circle), and Io (filled circle, on 
dashed circle) at western elongation are indicated. Io and the corotational plasma torus both 
move in the counter clockwise direction with the plasma torus overtaking Io with a relative 
period of ~13 hrs. The unmarked inner contour is 1.5 x 10 12 atoms cm 2 in (a) and 2.5 x 10" 
atoms cm 2 in (b). In (c), column density profiles are shown along the horizontal dashed line in 
(a) and (b). 

Figure 5. Pickup Ion Rates for Atomic Oxygen and Sulfur. The model calculated 
circumplanetary distribution of the pickup ion rate (ions cm 2 s" 1 ) integrated along the magnetic 
field is shown for electron impact ionization (El), charge exchange and (CE), respectively, in (a), 
and (b) as created by atomic oxygen in Figure 4a and in (d) and (e) as created by atomic sulfur 
in Figure 4b. The unmarked inner contour is 5 x 10 6 ions cm' 2 s' 1 in (a), (b), (d) and (e). In (c), 

pickup ion rate profiles for El and CE processes are shown along the horizontal dashed line in (a) 
and (b), respectively, while in (f), the same profiles are shown along the horizontal dashed line in 
(d) and (e), respectively. 

Figure 6. Net-Mass Loading Rate. The model calculated circumplanetary distribution for the 
net-mass loading of the plasma torus corresponding to the combined O and S pickup ion rates in 
Figure 5 is shown in (a) in two dimensions (kg m-2 s- 1 ), having been integrated along the 
magnetic field, and in (b) in one dimension (kg s-1 per 5°) as a function of System III longitude, 
having been integrated, in addition, along the radial coordinate. The relative locations of Jupiter 
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(half lit disk), Io’s orbit (dashed circle), and Io (filled circle, on the dashed circle) at western 
elongation are indicated. The two outer dashed lines in (a), originating from Jupiter, are centered 
20° on either side of the dashed line connecting Jupiter and Io. The locations of these same three 
dashed lines are also shown in (b) by the three vertical dashed lines and mark the boundaries 
between which 71% of the spatially-integrated net-mass loading rate of 154 kg s 1 is contained. 
In (a), the unmarked inner and outer contours are 3 x IO' 15 kg m 2 s' 1 and 2 x 10 18 kg m “ s ', 
respectively. 

Figure 7. Total-Mass Loading Rate. Same caption as Fig. Caption 6, except calculated for the 
total (El and equivalent CE) mass-loading rate, with the following changes. The inner contour in 
(a) is 8 x IO 15 kg m' 2 s' 1 . The ± 20° boundaries contain 71% of the total spatially-integrated 

mass-loading rate of 277 kg s' 1 . 

Figure 8. Mass per Unit Magnetic Flux Source Rate. The model calculated circumplanetary 
distribution for the mass per unit magnetic flux source rate (kg m ‘ s ’ T' 1 ) integrated along the 
magnetic field and corresponding to the combined O and S pickup ion rates in Figure 5 is shown 
in two dimensions in (a) and in a one-dimensional profile in (b) that is along the horizontal 
dashed line in (a). The unmarked inner contours near Io in (a) is 3 x 10 9 kg m' 2 s ’ T '. The total 

spatially integrated rate is 1.1 x 1 0 s kg s' 1 T 1 . 

Figure 9. Pickup Conductivity. The model calculated circumplanetary distribution for the 
height-integrated (integrated along the magnetic field) pickup conductivity (mho) for the 
combined O and S pickup ion rates in Figure 5 is shown in (a) and a one-dimensional profile in 
along the horizontal dashed line in (a) is shown in (b). The unmarked two inner contours near Io 
in (a) are 5 x 10 5 mho and 5 x 10' 4 mho. 

Figure 10. Radial Pickup Current. The model calculated circumplanetary distribution for the 
radial pickup conductivity (amps) for the combined O and S pickup ion rates in Figure 5 is 
shown in two dimensions in Jovicentric coordinates in (a) and in a one-dimensional profile in (b) 
along the horizontal dashed line in (a). 

Figure 11. Net-Energy Loading Rate. The model calculated circumplanetary distribution for 
the net-energy loading rate (watts m' 2 ) is shown in two dimensions in Jovicentric coordinates in 
(a) and in a one-dimensional profile in (b) along the horizontal dashed line in (a). The unmarked 
inner contour is 1 x 10 4 watts m' 2 . The spatially integrated net-energy rate is 1.03 x 10 1_ watts. 


37 



Atomic Oxygen Lifetimes 



4 5 6 7 8 9 10 11 

Plasmacentric Radial Coordinate (Rj) 


Figure 1 















Normalized Flux Distribution Function 



0 5 10 15 20 

Speed (kms'1) 


Figure 3 



Column Density (cm* 2 ) Jovicentric Perpendicular Distance (Rj) Jovicentric Perpendicular Distance (Rj) 



-10 -50 5 10 



Jovicentric East-West Distance (Rj) 


Figure 4 





Pickup Ions from Atomic Oxygen 


Pickup Ions from Atomic Sulfur 































Picki 



-10 -5 0 5 10 

Jovicentric East-West Distance (Rj) 


Figure 9 





Radial Pickup Current (amp) Jovicentric Perpendicular Distance (Rj) 



Figure 10 




Net-Energy Loading Rate (watts m-2) Jovicentric Perpendicular Distance (Rj) 



-10 -50 5 10 



Figure 11 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is est.mated to average 1 hour per response, including the time for reviewing 

and maintaining the data needed, and completing and reviewing the collection of information Send comments ragging this bu ^en e^.mate or ^ 

information, including suggestions for reducing this burden, to Washington Headquarters Services Directorate for formation ^ r ^ ns cmS?* ' 1215 JefferSOn D Hg ™ y ’ 
1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 


1 . AGENCY USE ONLY (Leave blank) 


2. REPORT DATE 
October 8, 2001 


3. REPORT TYPE AND DATES COVERED 
Final Report, Sept. 25, 1997 to Sept. 24, 2001 


TITLE AND SUBTITLE 

logenic Plasma and its Rotation-Driven Transport in Jupiter’s Magnetosphere 


5. FUNDING NUMBERS 


6. AUTHORS 


NASW-97023 


William H. Smyth 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Atmospheric and Environmental Research, Inc. 

131 Hartwell Avenue 
Lexington, MA 02421-3126 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

NASA Headquarters 
Headquarters Contract Division 
Washington, DC 20546 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


10. SPONSORING/MONITORING AGENCY 
REPORT NUMBER 



12a. DISTRIBUTION/AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

Model calculations are reported for the logenic plasma source created by atomic oxygen and sulfur above lo’s 
exobase in the corona and extended clouds (Outer Region). On a circumplanetary scale, two-dimensional 
distributions produced by integrating the proper three dimensional rate information for electron impact and 
charge exchange processes along the magnetic field lines are presented for the pickup ion rates, the net-mass 
and total-mass loading rates, the mass per unit magnetic flux rate, the pickup conductivity, the radial pickup 
current, and the net-energy loading rate for the plasma torus. All of the two-dimensional distributions are highly 
peaked at lo’s location and hence highly asymmetric about Jupiter. The logenic plasma source is also calculated 
on a much smaller near-lo scale to investigate the structure of the highly peak rates centered about lo s 
instantaneous location. The logenic plasma source for the Inner Region (pickup rates produced below los 
exobase) is, however, expected to be the dominant source near lo for the formation of the plasma torus ribbon 
and to be a comparable source, if not a larger contributor, to the energy budget of the plasma torus, so as to 
provide the necessary power to sustain the plasma torus radiative loss rate. 


14. SUBJECT TERMS 


15. NUMBER OF PAGES 
66 


Earth rotation, atmospheric angular momentum 


16. PRICE CODE 


17 SECURITY CLASSIFICATION OF 18 SECURITY CLASSIFICATION OF 19 SECURITY CLASSIFICATION OF 20. LIMITATION OF ABSTRACT 
REPORT THIS PAGE ABSTRACT 


Unclassified 


Unclassified 


Unclassified 


NSN 7540-01-280-5500 


Computer Generated 


STANDARD FORM 298 (Rev 2-89) 
Prescribed by ANSI Std 239-18 
298-102 










